Self sustained phantom wormholes in semi-classical gravity 

Remo GarattmQ 

Umversita degli Studi di Bergamo, Facoltd di Ingegneria, 
Viale Marconi 5, 24044 Dalmine (Bergamo) ITALY and 
INFN - sezione di Milano, Via Celoria 16, Milan, Italy. 

Francisco S. N. LobcQ 
Institute of Gravitation & Cosmology, University of Portsmouth, Portsmouth POl 2EG, UK and 
Centra de Astronomia e Astrofisica da Universidade de Lisboa, 
Campo Grande, Ed. C8 1749-016 Lisboa, Portugal 

A possible candidate for the late time accelerated expanding Universe is phantom energy, which 
possesses rather bizarre properties, such as the prediction of a Big Rip singularity and the violation 
of the null energy condition. The latter is a fundamental ingredient of traversable wormholes, and 
£ — ' it has been shown that phantom energy may indeed sustain these exotic geometries. Inspired by 

the evolving dark energy parameter crossing the phantom divide, we consider in this work a varying 
equation of state parameter dependent on the radial coordinate, i.e., w(r) = p(r)/p(r). We shall 
(N ■ impose that phantom energy is concentrated in the neighborhood of the throat, to ensure the flaring 

out condition, and several models are analyzed. We shall also consider the possibility that these 
phantom wormholes be sustained by their own quantum fluctuations. The energy density of the 
graviton one loop contribution to a classical energy in a phantom wormhole background and the finite 
one loop energy density are considered as a self-consistent source for these wormhole geometries. 
The latter semi-classical approach prohibits solutions with a constant equation of state parameter, 
which further motivates the imposition of a radial dependent parameter, w(r), and only permits 
\ solutions with a steep positive slope proportional to the radial derivative of the equation of state 

parameter, evaluated at the throat. The size of the wormhole throat as a function of the relevant 
■ parameters is also explored. 
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I. INTRODUCTION 



Unravelling the source for the late time accelerated expansion of the Universe has become a central theme in 
j^j), modern cosmology [if. Several candidates, responsible for this expansion, have been proposed in the literature, 
ranging from modified gravity theories to the introduction of cosmological models comprising of a negative pressure 
cosmic fluid, denoted as dark energy. The latter is parameterized by the equation of state to — p/p, where p is the 
spatially homogeneous pressure and p the dark energy density. A value of cj < — 1/3 is required to ensure the cosmic 
accelerated expansion, and lo = — 1 corresponds to the presence of a cosmological constant. It is interesting that recent 
constraints from observational data have confirmed that the dark energy parameter lies slightly below the characteristic 
cosmological constant value, to = — 1 @. Ameliorated fits to the data also indicate that an evolving equation of state 
parameter which crosses ui = — 1 is mildly favored [3j. The models satisfying lu < — 1 are denoted as phantom energy 
models, and possess rather interesting characteristics, in particular, the prediction of the Big Rip singularity and the 
violation of the null energy condition (NEC). The latter is a fundamental ingredient of traversable wormholes, and it 
was recently shown that phantom energy may indeed support these exotic geometries 0, H| ■ Inspired by the evolving 
dark energy parameter crossing the phantom divide, we consider in this work a varying equation of state parameter 
dependent on the radial coordinate, i.e., oj{t) = p(r)/p(r), and impose that phantom energy is concentrated in the 
neighborhood of the throat, in order to ensure the NEC violation, w(ro) < — 1, where ro is the wormhole throat radius. 
The latter restriction may be relaxed far from the throat. Despite the fact that the dark (phantom) energy equation 
of state represents a spatially homogeneous cosmic fluid, due to gravitational instabilities inhomogeneities may arise. 
Note that now one regards that the pressure in the equation of state p = ujp is a negative radial pressure, and the 
tangential pressure pt is determined through the Einstein field equations. This is fundamentally the analysis carried 
out in Ref. [6( , where the authors constructed a time-dependent solution describing a spherically symmetric wormhole 
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in a cosmological setting with a ghost scalar field. It was shown that the radial pressure is negative everywhere and 
far from the wormhole throat equals the transverse pressure, showing that the ghost scalar field behaves essentially as 
dark energy. Another interesting aspect of this topic is the effect of phantom energy accretion in wormholes immersed 
in this cosmic fluid. In fact, Gonzalez-Diaz analyzed the evolution of wormhole and ringhole spacetimes embedded in a 
background accelerating Universe [7[ driven by dark energy, and further considered the accretion of dark and phantom 
energy onto Morris-Thorne wormholes Q ■ It was shown that this accretion gradually increases the wormhole throat 
which eventually overtakes the accelerated expansion of the universe and becomes infinite at a time in the future 
before the big rip. 

Thus, the phantom wormhole models considered in this work may possibly arise from density fluctuations in the 
cosmological background @. In this context, we shall also be interested in the analysis that these phantom wormholes 
be sustained by their own quantum fluctuations. As noted above, wormholes violate the NEC, and consequently violate 
all of the other classical energy conditions. Thus, it seems that these exotic spacetimes arise naturally in the quantum 
regime, as a large number of quantum systems have been shown to violate the energy conditions, such as the Casimir 
effect. Indeed, various wormhole solutions in semi-classical gravity have been considered in the literature. For instance, 
semi-classical wormholes were found in the framework of the Frolov-Zelnikov approximation for (T^) [Ioj |. Analytical 
approximations of the stress-energy tensor of quantized fields in static and spherically symmetric wormhole spacetimes 
were also explored in Refs. [ill. However, the first self-consistent wormhole solution coupled to a quantum scalar 
field was obtained in Ref. [l 21 ] . The ground state of a massive scalar field with a non-conformal coupling on a short- 
throat flat-space wormhole background was computed in Ref. [13] , by using a zeta renormalization approach. The 
latter wormhole model, which was further used in the context of the Casimir effect )14J ]. was constructed by excising 
spherical regions from two identical copies of Minkowski spacetime, and finally surgically grafting the boundaries (A 
more realistic geometry was considered in Ref. [Hj]). In this context, this motivates the semi-classical analysis in this 
work, where the Einstein field equation takes the form 

G\iv = ft (Tjui/) ren 1 (1) 

with k = 87rG, and (T M „) ron is the renormalized expectation value of the stress-energy tensor operator of the quantized 
field. The formalism outlined in this work follows Ref. [l6j very closely. Now, the metric may be separated into a 
background component, g^ and a perturbation /i M „, i.e., g^ v — <? p „ + /i M „. The Einstein tensor may also be separated 
into a part describing the curvature due to the background geometry and that due to the perturbation, i.e., 

G l i V {g a p) = G^^ap) + AGfj,v(g a f3, h a p) , (2) 

where AG^^^q^, h a p) may be considered a perturbation series in terms of h^ v . Using the semi-classical Einstein field 
equation, in the absence of matter fields, one may define an effective stress-energy tensor for the quantum fluctuations 
as 

(T^y™ = -l(AG^(g a(3 ))™, (3) 

so that the equation governing quantum fluctuations behaves as a backreaction equation. The semi-classical procedure 
followed in this work relies heavily on the formalism outlined in Ref. [l6| . where the graviton one loop contribution 
to a classical energy in a traversable wormhole background was computed, through a variational approach with 
Gaussian trial wave functionals [l6l [l7| . A zeta function regularization is used to deal with the divergences, and 
a renormalization procedure is introduced, where the finite one loop is considered as a self-consistent source for 
traversable wormholes. Rather than reproduce the formalism, we shall refer the reader to Ref. [l6j for details, when 
necessary. In this paper, rather than integrate over the whole space as in Ref. [16], we shall work with the energy 
densities, which provides a more general working hypothesis. 

This paper is outlined in the following manner: In Section [TH we outline the field equations governing phantom 
wormholes with an r-dependent parameter, u)(r), and explore specific solutions. In Section HTT1 we analyze self- 
sustained traversable phantom wormholes in semi-classical gravity, considering that the equation governing quantum 
fluctuations behaves as a backreaction equation. Finally, in Section ITVl we conclude. 

II. FIELD EQUATIONS AND SOLUTIONS 

The spacetime metric representing a spherically symmetric and static wormhole is given by 

ds 2 = _ e 2*(,) dt 2 + dr* + ^ + ^2 e 
1 — Mr) r 
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where $(r) and b(r) are arbitrary functions of the radial coordinate, r, denoted as the redshift function, and the form 
function, respectively [l9j . The radial coordinate has a range that increases from a minimum value at r , corresponding 
to the wormhole throat, to infinity. A fundamental property of a wormhole is that a flaring out condition of the throat, 
given by (b — b'r)/b 2 > 0, is imposed [l9l [20j, and at the throat b(ro) — r — ro, the condition b'(ro) < 1 is imposed 
to have wormhole solutions. Another condition that needs to be satisfied is 1 — b(r)/r > 0. For the wormhole to be 
traversable, one must demand that there are no horizons present, which are identified as the surfaces with e 2 * — * 0, 
so that <i>(r) must be finite everywhere. 

Using the Einstein field equation, — with k — 8irG and c = 1, we obtain the following relationships 

b' = K r 2 p , (5) 

2r 2 (l-6/r) ' 1 ' 

Pr = -(Pt-Pr)-(p + Pr)&, (?) 

r 

where the prime denotes a derivative with respect to the radial coordinate, r. p(r) is the energy density, p r {r) is the 
radial pressure, and pt(r) is the lateral pressure measured in the orthogonal direction to the radial direction. 

Now consider the equation of state representing phantom energy, but with a varying parameter w(r), given by 
p r (r) = uj(r)p(r) with lu < — 1 in the phantom regime [H, 0]. The motivation for this analysis arises from a time- 
dependent crossing of the phantom divide, however, here we shall consider an r-dependent equation of state, with the 
imposition of phantom energy distribution in the throat neighborhood, with uj(ro) < —1, to ensure the NEC violation, 
and consequently the flaring-out condition. This condition may be relaxed far from the throat. Taking into account 
Eqs. ([H])-©, we have the following condition 

b(r)+u(r)rb'(r) 

( > 2r 2 [l-6(r)/r] ' 1 ' 

We now have four equations, namely the field equations, i.e., Eqs. ©-([T]), and Eq. ([5]). with six unknown functions 
of r, i.e., p(r), p r {r), Pt(f), b(r), <&(r) and w(r). One may adopt different strategies to construct solutions with the 
properties and characteristics of wormholes. One may consider a plausible matter-field distribution, by considering 
convenient stress-energy tensor components, and through the field equations and Eq. ([5]) determine the metric fields 
and the function w(r). 

Alternatively, by appropriately choosing the form function and ui(r), one may integrate Eq. ((SJ), to determine the 
redshift function, $(r). As an example, consider the following functions 



b(r) = r + 7^r (l-f) , (9) 
oj( r ) = -a 2 + p 2 (l - ^) , (10) 

with < 7 2 < 1, so that one obtains a wormhole solution with b'(r) = ^r^fr 2 , and at the throat we have b(ro) = ro 
and b'(ro) < 1. Relatively to the function cu(r), we impose u>(ro) = —a 2 < — 1, and —1 < u>(r) = —a 2 + f3 2 < —1/3 
as r — * oo. Now, substituting Eqs. (f9j)- (fT0]) into Eq. J5J), we obtain the solution 

*(r) = - W ^ ) Mr-r Q )--(l + a 2 + -; 



This corresponds to a non-traversable wormhole solution, due to the presence of an event horizon, as may be readily 
verified from the first term in the right hand side. Therefore, we fine-tune the parameter a 2 = I/7 2 , so that the 
redshift function finally takes the form 

which is now finite throughout the range of r, and where we have considered C = 0, without a significant loss of 
generality. Note that this solution corresponds to an asymptotically flat traversable wormhole, as b(r)/r — > and 
$(r) — > for r — > 00. 
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One may also consider restricted choices for $(r) and oj(r), and through Eq. 
through the following relationship 



b(r) = e 



_ o -T(r ,r) 



ro 



with 



r(r ,r) = 



2$'(f)f + l 



determine the form function 



(13) 



uj{r)r 



dr . 



(14) 



Now, interestingly enough results may be inferred by imposing a constant redshift function, $' = 0, i.e., b(r) 
Lo(r)rb'(r) = 0, which may be integrated to yield the following form function 



b(r) = tq exp 



dr 



uj{r)i 



(15) 



The latter may also be deduced from Eq. (|13|). by setting $'(r) = 0. Note that for a constant equation of state 
parameter, lo, the form function is given by b(r) = ro(r /ro) _1/ '' a, , which is a solution analyzed in detail in Rcf. 
For instance, consider the following function 



w(r) 



a 2 (3 2 (r/r ) 



a 2 (1 - r/r ) - /3 2 



which takes the values w(ro) = —a 2 < —1 at the throat, and cu(r) — —j3 2 as r 
function in Eq. (|15|) . we deduce the following form function 



(16) 



oo. Now, substituting this specific 



b(r) 



r 

ro | — 

ro 



1//3 2 



exp 



(17) 



Note that if we wish to impose the crossing of the phantom divide the condition 1 < 1 needs to be satisfied. However, 
from the expression of the form function, we verify that in order to have an asymptotically flat spacetime, i..e, b/r — > 0, 
as r — ► oo , the condition f3 2 > 1 is imposed. Therefore for this particular choice of w(r), the solution remains in the 
phantom regime for all values of r, so that the wormhole is constituted totally of phantom energy. The derivative 
of the form function evaluated at the throat, is given by b'(ro) = 1/a 2 < 1, consequently satisfying the flaring-out 
condition. 

Now, consider the function given by Eq. (|10p . then integrating Eq. (|15p. we deduce the following form function 



b(r) 



ro 



r/3 2 



r 



1 - 



P 2 



l/(« 2 -/3 2 ) 



(18) 



Note that for this case to be an asymptotically flat spacetime, then we need to impose a 2 — 1 > 1. This latter 
restriction imposes that u> remains in the phantom regime, w(r) < —1 for all r. The derivative of the form function, 
evaluated at the throat, is given by b'(ro) = 1/a 2 < 1, which also satisfies the flaring out condition. We note that it 
is possible, in principle, to construct a wide variety of traversable wormholes supported by phantom energy, with an 
r-dependent equation of state parameter. 



III. SELF-SUSTAINED PHANTOM WORMHOLES 

In this Section we shall consider the formalism outlined in detail in Ref . [l6[ , where the graviton one loop contribu- 
tion to a classical energy in a wormhole background is used. The latter contribution is evaluated through a variational 
approach with Gaussian trial wave functionals, and the divergences are treated with a zeta function regularization. 
Using a renormalization procedure, the finite one loop energy was considered a self-consistent source for a traversable 
wormhole. We refer the reader to Ref. [IB] for details. 

The classical energy is given by 
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where the background field super-hamiltonian, 7i^°', is integrated on a constant time hypersurface. R is the curvature 
scalar, and using metric (|1]), is given by 



R 



-2| 1-- 

r 



$" + ($ 



/\2 



36 -Ar 



r[r — b) 2r(r — b) 



We shall henceforth consider a constant redshift function, $'(r) = 0, which provides interestingly enough results, so 
that the curvature scalar reduces to R = 2b' /r 2 . Thus, the classical energy reduces to 



2G 



dr r 



b(r) 



^Jl-b{r)/r uj{r)r 3 ' 



(19) 



where the condition b'(r) = —b(r)/[ui(r)r] has been used to eliminate the b'(r) dependence. 
We shall also take into account the total regularized one loop energy given by 



E TT = 2 



dr- 



y/1 -b(r)/r 



[pi{e) + p 2 (e)] , 



(20) 



where once again, we refer the reader to Ref. [16] for details. The energy densities, Pi(e) (with i = 1, 2), are defined 
as 



4?r Jy/lW) 



dE, 



E— 1/2 



64tt 2 



2 In 2 - - 
2 



(21) 



The zeta function regularization method has been used to determine the energy densities, pt. It is interesting to 
note that this method is identical to the subtraction procedure of the Casimir energy computation, where the zero 
point energy in different backgrounds with the same asymptotic properties is involved. In this context, the additional 
mass parameter p has been introduced to restore the correct dimension for the regularized quantities. Note that this 
arbitrary mass scale appears in any regularization scheme. 

Note that in analyzing the Einstein field equations, one should consider the whole system of equations, which in 
the static spherically symmetric case also includes the rr-component. The joint analysis of the equations is necessary 
to guarantee the compatibility of the system. However, in the semi-classical framework considered in this work there 
is no dynamical equation for the pressure. Nevertheless, one may argue that the semi-classical part of the pressure is 
known through the equation of state that determines the relation between the energy density and the pressure. 

Since a self sustained wormhole must satisfy 



H^ = -E TT , 



(22) 



which is an integral relation, this has to be true also for the integrand, namely the energy density. Therefore, we set 

1 b(r) 



2Guj{ry 



= -2 [ Pl (e) + p 2 (e)} 



(23) 



At first sight one can think that the phantom region is straightforward to analyze. However, within the approach of 
Ref. [la ], this is not trivial as the main difference is in the integration. In this paper, we shall consider the approach of 
working with the energy density, while in Ref. (To) , the integration over the whole space was considered. In Appendix 
IA1 we show that the classical term integrated over all space forbids the appearance of phantom energy [l8| . at least 
for constant uj. 

For this purpose, the Lichnerowicz equations provide the potentials, which are given by 



Ui(r) 



G 



1 



Hr) 

r 
bjr) 
r 



3 

2^2" 

36(r) 
2r 3 



b'(r) 



u(r) 



6(r) 
r 

1 



(24) 



U 2 (r) 



v 
G 



6(r) 

r 
6(r) 

r 



-b'(r) b(r) 



3 

2^ 
36(r) I" 1 
2r 3 |_3w(r) 



(25) 
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where the relationship 6'(r) = —b(r)/[uo(r)r] has been used once again in the last terms, to eliminate the b'(r) 
dependence. We refer the reader to Ref. for the deduction of these expressions. 
Thus, taking into account Eq. (|2"Tjl . then Eq. (|23|) yields the following relationship 



1 b(r) 



1 



2G uj(r)r 3 32n 2 s 



[Ul{r)+U 2 2 {r)} 



32tt 2 





4/i 2 




4/i 2 


)] 








u 2v ^ 





(26) 



It is essential to renormalize the divergent energy by absorbing the singularity in the classical quantity, by redefining 
the bare classical constant G as 



1 



2 [UUr) + Ui(r)) c(r)r 3 



Using this, Eq. (|26]) takes the form 

1 b(r) 



2G w(r)r 3 



32tt 2 





32tt 2 


b(r) 








4 M 2 




4/i 2 


)] 













(27) 



(28) 



Note that this quantity depends on an arbitrary mass scale. Thus, using the renormalization group equation to 
eliminate this dependence, we impose that 



dfi 



1 



b(r) 



2G uj(r)r 3 



fj, 



32tt 2 dn 





4/i 2 


) + U% In ( 


4/x 2 


)] 








c/ 2 ^ 





which reduces to 



(29) 



(30) 



The renormalized constant Go is treated as a running constant, in the sense that it varies provided that the scale /i 
is varying, so that one may consider the following definition 



1 



G (/i) G (/i ) 
Thus, Eq. (|26|) finally provides us with 



1 , [Ul(r) + Ui(r)} w (r)r 3 ^ ( p 



1 



ui(r)i 



G (/i ) 16^ 2 6(r) 



£/ 2 (r)ln 



8tt 2 



£/i(r)VS 



6(r 
+ C/ 2 (r)ln 



Mo 



U 2 {r)^~e 



(31) 



(32) 



Now, the procedure that we shall follow is to find the extremum of the right hand side of Eq. ([321) with respect to r, 
and finally evaluate at the throat ro, in order to have only one solution (see discussion in Ref. [161]). 

To this effect, we shall use the derivative of the potentials, Eqs. J24j) - (l2"5|) . which we write down for self-completeness 
and self-consistency, and are given by 
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Ui(r) = --t 1-- 



r 
12 



6(6 - b'r) 3 

r 4 + 2 

6(6 - b'r) 3 

7* + 2 



36 /l 



- + 1 I --^T 

U) 



b'r -36/1 



3lJ 



b u' 
r 3 3cj 2 



and once evaluated at the throat take the following form 



'o 
3K 



with 
with 



A = 2 
B = 2 





1 


^0 / 




l + w \ 


1 


Wo / 





[(l + 3w )(l+w )+o; r ] 



(1 + 3cj )(1-3w ) + 



w r 



(33) 
(34) 

(35) 
(36) 



where the relationship 6' = —b/(ujr) has been used to eliminate the 6' dependence. We have defined w(ro) = uiq, and 
the term u>' ro is a number as u>' oc 1/ro. The potentials, Eqs. (|24|) - ([25|) evaluated at the throat reduce to 



3Y_ 

2 
(] 



with X = 
with Y = 



1 + ^0 

3lu — 1 
3u) 



(37) 
(38) 
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To determine the extremum of Eq. (j32[) with respect to r, we shall consider, for convenience, the following definitions 



tu(r)r 3 



F{r) = _[i + 3w(r)+w'(r)r], 
b(r) 



(39) 



where in the latter expression, b' = —b/(ur) is used. Thus, the extremum of Eq. (|32p with respect to r, takes the 
form 



Hr) 



2UxU' x hx 



Ui\e 
4M§ 



2U 2 U!> In 



gaje 



0. 



where the following relationship 



Win 



4^ 



(40) 



(41) 



has been used. 

Now, evaluated at the throat, ro, in order to have only one solution, using Eqs. (|35p - (l36|) . and the expressions for 
the potentials, Eqs. (|3"7] ) -(|58 )) . we verify that Eq. (@DJ) provides 



-(1 + 3^o r o) 



Wlrl 



- w 



XA In 



/ 3eX 
Wort 



YB ln 



/ 3eY 
Wor 2 o 



0. 



Going through the whole tedious, but fairly straightforward, calculation of the throat, we arrive at 



8nl 



where the constants are defined as 



A+C B + D I A + B 

il , with U — X e Y E exp 

2E 



A = uj XA , B = -uj YB , 
C= i(l + 3w +4r )X 2 , 

D = + 3w +w r )F 2 , 



(42) 

(43) 
(44) 

(45) 



and i? = A + B + C + D. We may now consider uio < — 1 (phantom regime) and w ro as free parameters, and 
consequently plot ro , which is depicted in Fig. Q] We have defined the dimensionless quantity Rq = fo/ (^J §^?j = V 7 ^, 
which is depicted by the surface in Fig. [T] Note that the throat size decreases for increasing values of the parameter 
X = uj' Q r . 




20 -1 



FIG. 1: We have defined the dimensionless quantity Ro = ro/ ( \/§^r ) = which is depicted by the surface. Note that the 



8n 

throat size decreases for increasing values of the parameter x = Wo^o- 



Now, restrictions on the parameter range of lvq and w ro may also be imposed. The quantity Gq(jj,o) must be 
positive, so that Eq. (|3"2")l evaluated at the throat, reduces to 



1 



^ 3|^o| 



ln 



X\ X (Y- Y 



ft ) 



(46) 
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The log must be positive, which places restrictions on the range of iv' r . We shall depict F = ln[(X/n) x (Y/n) Y } 
in Fig. [5] Note that only solutions with a steep positive slope proportional to lu' , at the throat, are permitted. 
As loq — > —1, we verify that w ro increases significantly. Solutions with a constant parameter u, i.e., lu' q — 0, are 
prohibited, which is in agreement with the calculations in Appendix |A"1 




FIG. 2: The log in the definition of Go(/io) must be positive, which places restrictions on the range of x = u)' ro. The surface 
in the figure is given by F(x,u>o) = hi[(X/Q) x (Y/Q) Y ]. Only solutions with a steep positive slope of the parameter uj'qTo, at 
the throat, are permitted. As too tends to — 1, we verify that u}' ro increases significantly. Note that solutions with a constant 
parameter w, i.e., u' = 0, are prohibited. 



An interesting generic constraint on [ij \i§ may also be deduced. Substituting Eqs. 
evaluated at the throat, provides the following relationship 



1 



= 3|^ |mq 



In 



X\ X (Y^ Y 



9. 



2(X 2 +Y Z ) 



and into Eq. (3TJ, 
(47) 



Now Go (a*) must be positive, so that the above relationship places a constraint on /i//io, namely 



A* < Mo 



n 



n 



(48) 



The above results have been based on a variational calculation, concerning the one-loop graviton contribution in 
a wormhole background. On general grounds, in order for the one-loop contribution to stabilize the wormhole it 
can only do so at the Planck scale, so that the wormhole throat is expected to be of the Planck order. We have 
found an upper bound on the mass parameter fx, and we emphasize that we are not interested in computing an exact 
result of the latter, but in the possibility that phantom wormholes be, in principle, sustained by their own quantum 
fluctuations. Indeed, we have demonstrated this fact in the context of an r-dependent equation of state parameter 
within the one-loop graviton contribution approach. 



IV. CONCLUSION 



Inspired by the evolving dark energy parameter crossing the phantom divide, we have considered in this work a 
varying equation of state parameter dependent on the radial coordinate, i.e., u>(r) — p(r)/p(r), and imposed that 
phantom energy is concentrated in the neighborhood of the throat. The pressure in the equation of state is a negative 
radial pressure, and the tangential pressure pt is determined through the Einstein field equations. We have also 
considered the possibility that these phantom wormholes be sustained by their own quantum fluctuations. The 
energy density of the graviton one loop contribution to a classical energy in a phantom wormhole background, and 
the finite one loop energy was considered as a self-consistent source for these phantom wormholes. The size of the 
wormhole throat as a function of the parameters luq and w r o was further explored, and in order for the one-loop 
contribution to stabilize the wormhole it can only do so at the Planck scale, so that the wormhole throat is expected 
to be of the Planck order. In a rather speculative scenario, one may also argue that the accretion of phantom^nergy 
by the wormhole gradually increases the wormhole throat to macroscopic size, much in the spirit of Refs. 0, [H- 

At first sight one can think that the phantom region is straightforward to analyze. However, within the approach of 
Ref. [161 ], this is not trivial. Indeed, the main difference is in the integration. In this paper, we have worked with the 
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energy density, while in Ref. [T(|, the integration over the whole space was considered. In Appendix Kl we show that 
the classical term integrated over all space forbids the appearance of phantom energy for constant u> [18| . However, it 
is important to emphasize that the semi-classical approach used in this work only rules out self sustained wormholes, 
with constant u>, concerning the one-loop graviton contribution, and not the whole plethora of phantom wormholes, 
as the ones analyzed in Refs. [HQ. A few words should be spent concerning the difference between the phantom and 
the dark energy case. Actually the dark energy (not energy density) could be considered, because the classical term 
strongly forbids phantom energy. Nevertheless, the dark energy region is completely incompatible with asymptotic 
flatness [ijj]. The analysis for co(r) is not a trivial task and it depends on a case to case. However, we have shown that 
interesting solutions do, in fact, exist, namely those that permit solutions with a steep positive slope proportional to 
the radial derivative of the equation of state parameter, evaluated at the throat. It is rather important to emphasize 
a shortcoming in the analysis carried in this paper, mainly due to the difficulties of the semi-classical theory of gravity 
and, in particular to the technical problems encountered. We have considered a semi-classical variational approach 
using the energy density, which imposes a local analysis to the problem. Thus, we have restricted our attention to 
the behavior of the metric function b(r), through a variational approach, at the wormhole throat, ro. Despite the 
fact that the behavior is unknown far from the throat, due to the high curvature effects at or near ro, the analysis 
carried out in this work should extend to the immediate neighborhood of the wormhole throat. Therefore, we may 
conclude that, in principle, phantom wormhole throats with a varying equation of state parameter dependent on the 
radial coordinate, may be self sustained by their quantum fluctuations within the semi-classical approach outlined in 
this paper. 
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APPENDIX A: COMPUTATION OF THE CLASSICAL TERM 



Here, we compute the hamiltonian of Eq. (|19[) when u>(r) is a constant. It is easy to show that 



i 

In the previous integral there is an extra factor "2" coming from the counting of the universes. We change variable 
to obtain 

2r f +co dt _ r (1 1 \ _ r , 

Gcox J C0S h 2 - 2 /- (t) ~ G ( 1 + co) B \ 2 ' 1 + co J- G A M ( A2) 

with 







(l+w) 




R r 




(1- 


f w)r( 


' 3+lo \ 




^2+2u> J 



In Eq. (|A2[) , we have used the following formula 

cosh" t 2 \ 2 ' 2 J \Re(^-^)<0 v ' 

The existence of the integral is submitted to the following condition (/i = 0) lo > — 1 and this shows that the phantom 
region is forbidden. 
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